Dispersionless motion in a driven periodic potential 
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Recently, dispersionless (coherent) motion of (noninteracting) massive Brownian particles, at 
intermediate time scales, was reported in a sinusoidal potential with a constant tilt. The coherent 
motion persists for a finite length of time before the motion becomes diffusive. We show that such 
coherent motion can be obtained repeatedly by applying an external zero-mean square-wave drive 
of appropriate period and amplitude, instead of a constant tilt. Thus, the cumulative duration 
of coherent motion of particles is prolonged. Moreover, by taking an appropriate combination of 
periods of the external field, one can postpone the beginning of the coherent motion and can even 
have coherent motion at a lower value of position dispersion than in the constant tilt case. 

PACS numbers: : 05.10.Gg, 05.40.-a, 05.40. Jc, 05.60.Cd 



The inertial Brownian particle motion in periodic 
potentials [l|,l2| has been an archetypal model to theoret- 
ically understand many phenomena in physical systems. 
The current-voltage characteristics of (RCSJ model of) 
Josephson junctions[3], the electrical conductivity of su- 
perionic solids motion of adatoms on the surface of a 
crystal etc., are some of the important exam pies [J. 
However, not all behaviour of the model particle mo- 
tion in all time regimes are exhaustively investigated. 
A recent example being the discovery of dispersionless 
particle motion in a tilted periodic potential in the in- 
termediate time regime by Lindenberg and coworkers @. 
During the coherent motion, the ensemble averaged posi- 
tion dispersion, Ax(t) —< (x(t)— < x(t) >) 2 >, remains 
constant. 

This interesting phenomenon is shown (numerically) 
by particles moving on a cosinusoidal potential with a 
constant tilt (CT), i<o, in a medium with constant fric- 
tion coefficient 70 [6| in a limited (-Fb>7o) region. The 
particles, after crossing the immediate barrier move (after 
t = T\ > tki the Kramers mean passage time) coherently 
with velocity v ~ — . The coherent motion continues un- 
til it is overwhelmed (at around t — T2) by the diffusive 
motion of the particles. t\ and T2 are specified only as a 
rough guide 6]. In this work we investigate the effect of 
a zero-mean square- wave external drive (ZMSW) F(t) of 
half period r and amplitude Fq, instead of a CT. 

The coherent motion, naturally, gets interrupted upon 
reversal of direction of F at t — t (t\ < t < T2). Inter- 
estingly, as a main result of this work, the coherent mo- 
tion once disturbed, by reversing the field at t = r, gets 
reestablished around t — r + t± in almost the same form 
as it was during t± < t < t-i in the CT case. And this 
loss and subsequent recovery of coherent motion contin- 
ues for a large number of reversals of F(t) . The dispersion 
Ax(t), however, increases rapidly during r < t < t + ti. 
(During < t < n, Ax(t) ~ i Q , a » 2). 



We consider the motion of an ensemble of Brow- 
nian particles each of mass m moving in a poten- 
tial V(x) = —Vq sin(fca;) in a medium with friction 
coefficient [7| j(x) = 7o(l — Xsm(kx + </>)) at temper- 
ature T (in units of ks) and subjected to an external 
force field F(t). The corresponding Langevin equation in 
dimensionless form is given byQ 

^ = -l(x)—+cosx + F(t) + ^(xyTat), (0.1) 

with 7(x) = 7o(l — Asin(a; + <fi)). Here all terms are 
written in units of m, Vq and k, for example, the scaled 
time t is measured in units of 



etc. The Gaus- 
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sian distributed zero-mean random forces £(t) satisfy 
< £(t)£(t') >= 25 (t — t'). The applied square-wave force 
F(t) is taken as F{t) = ±F for (2nr < t < (2n + l)r) 
and F(t) = ^F for ((2n + 1)t < t < (2n + 2)r) with 
n = 0,1,2,.... For constant applied force F(t) = Fq 
(for all t) the equation is solved using the matrix con- 
tinued fraction method and also numerically, support- 
ing each other quantitatively 0, d, [TH EH E3- However, 
for finite r, the equation could be solved only numer- 
ically. The integration of the equation, using the 4 
order Runge-Kutta method [l3j], was carried out in time 
steps of At = .001. 

We take A = 0.9. A 7^ 0, however, is relevant only while 
discussing ratchet current at the end. In fact, except for 
this minor point, all the results discussed in the following 
qualitatively remain same for the simpler case of A = 0. 

The particles exhibit coherent motion in the potential 
V(x) — — sm(x) — xFo, in the intermediate times roughly 
in the range [ti(« 2 x 10 3 ) < t < r 2 (« 3 x 10 4 )], for 
F = 0.2. Therefore, we choose the ZMSW field F(t) of 
amplitude Fq = 0.2 and, in most cases, half period r = 
5000 which is well within the range [ti,T2]. Naturally, in 
the first half period (0 < t < t) the motion is same as 
in the CT case. In all cases, we take the initial (t = 0) 
particle position distribution as 5(x — §) and Maxwell 
velocity distribution corresponding to T — 0.4. Note 
that after every half period r the periodic potential gets 
tilted in the reversed direction as a result of field reversal. 
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FIG. 1: The position dispersion Ax(t) for r = 5000 (a), 400 
(b), and 200 (c) are plotted. The Inset shows curve (a) ex- 
tended to 10 half periods. 
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FIG. 2: Ax(t) for three r = 5000 (a) ten r = 400 and two 
r = 5000 (b), five r = 200 and two r = 5000 (c), five r = 200 
and one r = 10000 (d), and one r = 5000 and five r = 400 
(e) half periods are plotted. 



Fig.l, curve (a) and its extended plot in the Inset show 
that by applying the field, F(t), with r = 5000 a repet- 
itive sequence of trains of coherent motion, with char- 
acteristic constant Ax(t) is obtained. These bursts of 
coherent motion are quite robust. Each burst of coher- 
ent motion is preceded by a length of dispersive particle 
motion. As a result Ax(t) grows, in discrete steps, with 
time, as the number n of half periods increases. 

The generation of coherent motion continues for many 
(n >> 8) half periods (r = 5000) of F(t), Inset of Fig.l. 
Thus, as an important consequence, the cumulative du- 
ration (» n(r - ti)) of coherent motion when driven by 
F{t) is made much larger than the duration, 72 — n, 
achievable in the CT case. 

In the Fig.l are also plotted Ax corresponding to the 
half periods r =400 (curve (b)), and 200 (curve (c)). The 
important feature to be noticed in the figure is that, for 
t = 400,200 < tk,ti, Ax dips immediately after the 
field is reversed before it rises. This behaviour of dis- 
persion dipping and subsequent rise becomes most pro- 
nounced at a small but intermediate r. It continues for 
many periods of F(t). 

Fig. 2, summarizes the main results of this work. It 
shows many interesting effects of taking few initial half 
periods of F(t) of smaller duration, r = 400 (curve (b)), 
and r = 200 (curves (c-d)) and then making its later 
half periods r = 5000 or larger (e.g., 10000, curve (d)). 
The curves (b-c) show that, in this case too, during the 
later half periods r = 5000 of F(t), a similar sequence 
of trains of coherent motion, as in Fig.l, Inset (or curve 

(a) ), can be obtained. It also allows to postpone (curve 

(b) ) the appearance of coherent motion beyond t = t\. 
Moreover, it is possible to obtain coherent motion with 
lower constant dispersion (curves (c,d)) than in the CT 
case (i.e., lower than the constant Ace in the first half 
period in the curve (a)) too. Curve (e) shows that the 



curves (b-d) can be repeated, using the same customized 
procedure, many times over again. 

The results of Figs.l and 2 can be understood by 
analyzing the time evolution of velocity distribution, 
P(v). P(v) at various phases of F(t) at t = 5r, and at 
t = t + 15.6 for r = 5000 and at t = 2t, and t = 2t + 16 
for r = 1000 are shown in Fig. 3. P(v) invariably assumes 
almost a Gaussian form of same width and centred at a 

fixed v i=a ±— at t — nr for t = 5000. The Inset, show- 

70 

ing the mean velocity v(t) and velocity dispersion Av(t), 

supports this observation. However, at t = nx 1000, P(v) 

has two peaks one centred at v « ±— and the other at 

~ 70 
v = 0. 

The P(v) peak at v = shows that at t — nr, for 
r < ti , some particles are left behind in the locked state 
in one or some other wells. The smaller is t more promi- 
nent this latter peak is left at t = t. These locked par- 
ticles try to gravitate to their respective well bottoms 
and, being slower, even succeed in shrinking the position 
distribution P(x) more effectively than the much faster 
particles in the running state on the front. Thus, for 
smaller t < ti, the bimodal nature of P(v) helps appre- 
ciable Ax dipping immediately after field reversal, Fig.l. 

At t = nr + 1 , 15 < to < 16, v becomes zero for all t, 
and n. (It turns out, 2to ~ ^.) The corresponding P(v) 
are also shown in Fig. 3. For t = 5000, P(v) at t = nr+t 
is like a sum of two "Gaussians" centred on either side 
of v = 0. But for t = 1000, the surviving P(v) peak at 
v = at t = nr, contributes an additional one centred at 
v = 0, making v = at the same delay time To- 

Due to the tilt direction reversal the particles, in the 
running state, are forced to reverse their direction of mo- 
tion and hence each of them necessarily go through zero 
velocity at least once momentarily. Thus, the entire sys- 
tem passes approximately through a "thermal" state at 
t = nr+To- Hence, after every field reversal at t = nr the 
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FIG. 3: P(v) at t = 25000 (right peak) and at 25015.6 (middle 
thin line) for r = 5000, and at t = 2000 (bimodal), at 2016 
(middle bold line) for r = 1000. v(t), and Av{t) for r = 5000 
are plotted (v = ±-fo/7o lines drawn) in the Inset. 
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FIG. 4: Ax(t) for r = 1000 (a), 5000 (b), and for the case 
of CT Fo = 0.2 (c) averaged over, respectively, 20, 60, and 
18 ensembles are plotted. Lines of slope 1 are fitted to the 
curves. The short line at the lower left corner indicates Aa; ~ 
t a ,a m 2. The inset shows variation of D with r. 



particles begin their subsequent journey in the reversed 
direction with almost the same delayed initial condition 
(at t — tit + tq) of thermalized P(v). Therefore, Aa; are 
expected to behave similarly after every nr. 

As discussed above, it is just the reversal of field which 
leads P(v) to the required form at t = nr+to, irrespective 
of the form of P(v) at t = nr for any value of r >> To ~ 
16. It shows that in order to obtain coherent motion 
neither an exact initial Gaussian velocity distribution is 
necessary nor all the particles are required to be intially 
confined sharply to a single well bottom of the periodic 
potential. Also, a mere switching the field alternately on 
(Fq 7^ 0) for duration r and off (Fo = 0) for the same 
duration r, fails to yield results like Figs.l and 2. This 
indicates that a reversal of the field direction is essential 
because this alone ensures a "thermalized" P(v). 

It must also be noted that during the CT case the 
average particle displacement is large, w by the end 
of its coherent motion whereas in the ZMSW case it is 
zero, for A — 0, and small and finite for A / and 
4> 7^ 0, 7r after any large time t = 2nr. This is an added 
practical advantage over the CT case for, in the ZMSW 
case, most of the particles on the average remain confined 
to a finite region of space despite periodically showing 
coherence of motion for a long time. 

The dimensionless values of parameters used (e.g., 
7o = .035, T = A,t = 5000, and F = .2) and other 
derived quantities when restored to their usual units are 
presented in Table 1 for three illustrative cases: (i) The 
motion of an Ag ion in Agl crystal[14], (ii) the mo- 
tion of a macromolecule (kinesin) along a polymer fi- 
bre (microtubule) [IH] , and (iii) diffusion of Cooper pairs 
across a Josephson junctional. Notice that cu(— << 
u>o, in the particle motion case, and u << up, the 
Josephson plasma frequency (= showing that the 



systems considered are, indeed, underdamped. The last 
column of the Table gives the magnitude of mean velocity 
(mean voltage) attained during the coherent state when 
the initial value of the drive field F(t) (I(t)) is fixed ei- 
ther at + \F 1 (+1-^0 1) or with their sign reversed and not 
an equal mixture of both. Also, during the half period r 
the particles move to an average distance (the product of 
quanties in the last two columns) of the order of a micron 
(fi) which will get retraced in the next half period. This 
gives a rough idea of the sample size one would need to 
take in a ZMSW case. Also, ^ << w (up) (by about 
two orders of magnitude). 

The calculated average velocities v and velocity dis- 
persions Av are plotted, in the Inset of Fig. 3. An equal 
mixture of F(t = 0) = ±|-Fo| makes v close to zero during 
coherent motion. It is exactly zero for A = at all t and 
hence even in the limit t — > 00 v remains zero. However, 
for A 7^ and <p 7^ 0, 7r a nonzero finite mean (steady 
state) velocity is obtained[l(| HH, 03 earlier. The con- 
tribution of coherent particle motion being insignificant, 
the dispersive motion alone contributes to this ratchet 
current of particles. 

The diffusion constant, D, defined as liiut_ ) . 00 Ax(t) = 
2Dt, is hard to calculate for a constant tilt Fop. The 
asymptotic limit barely reaches even by t = 10 7 . How- 
ever, for ZMSW this limit is readily reached by t = 10 7 , 
Fig. 4. It may be noted that for each curve in Figs.l and 2 
we have averaged over 2000 realizations but in Fig. 4 we 
could average over number of realizations ranging only 
between 18 and 60. The nature of Ax(t) so clear in Fig.l 
appears less convincing in Fig. 4. Therefore, it is hard to 
conclude that the same nature of Ax, as in Fig.l, will 
continue till the asymptotic time regime. However, from 
the "thermalized" P(v) argument given earlier, there is 
a fair likelyhood that the nature of Ax(t) shown in Fig.l 
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TABLE I: The systems considered are: (1) Ag: Ag on Agl lattice, (2) Mm: a macromolecule along a polymer, and (3) JJ: 

Josephson junction. The symbols have their usual meaning. The RCSJ-model JJ equation equivalent to Eq.(2.1) is[3|: 

|f + £G(1 + Acos(6> + (/>))§ + h sin(0) = I(t) + ^/2TG(l + Acos(0 + <j>))£(t). (Add -§ to 6 for exact correspondence.) 
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will extend to a large number of half periods of F(t) 
provided a large number of particles are considered for 
averaging. 

The nature of Ax(t) shown by the curves in Figs.l, 2 is 
in no finite region close to Ax(t) ~ t. It remains, there- 
fore, open to explain why a large number of repeatedly 
same diverse combinations of dispersions such as ones 
ranging from Ax(t) ~ t 2 to Ax(t) ~ i°, when averaged 
over a large number of realizations, yields the same na- 
ture of dispersion, lim t ^ocAx(t) ~ t, for all t, Fig. 4. 

D(t), plotted in the Inset of Fig. 4, are rough estimates 
as the averagings are done only over a small number of 
realizations. However, the overall qualitative trend of 
D(t) remains valid. D(t) shows a peak around r = t\. 
For r > t\ , the closer t is to t\ smaller is the constant Ax 
region and hence larger is the fraction of sharply rising 
Aa; region. Naturally total Ax will be larger as r — ► t\. 
However, the nature of D(t) as r — * T2 is not clear from 
the available data. In the range r < Ti, Ax rises only 
after an appreciable dipping, Fig.l. Therefore, the initial 
rise of Aa; is slower as r is decreased from t\ resulting in 
a smaller total Ax(t) as t — > oo and hence smaller D. 

From the rms spread (V Ax) point of view the advan- 
tage of ZMSW F(t) over CT, except in cases like curve 
(d) in Fig. 2, quickly evaporates as t increases. Whereas 
(for an Ag particle) at t = lOr = 12.3 x lCT 9 s for the CT 
(ZMSW) case, the mean displacement x is 1.8/x (ss 0) and 



VAx = .14/i (.40//), at t = 2.46 x 10 6 s the correspond- 
ing numbers are x — 3.64mm (w 0) and V Aa; = .77/i 
(6.08/i). Perhaps in the ZMSW case, the particles left 
behind during a r get pushed father away during the 
next r, make the \J Ax(t) increase faster as t increases. 

Coherent motion is observed only in the negative slope 
region of D(r). However, for this same system it is shown 
in Ref.[l2j that ratchet current is maximum for a value 
of r ~ 500, i.e., in the rising D(t) region and becomes 
significantly small for larger r > t% and almost zero at 
r = 5000. The peak of the D(t), thus, roughly divides r 
into two regions: (i) small r giving ratchet current, and 
(ii) larger r showing coherent motion. 

To summarise, the dispersionless (coherent) motion 
discovered earlier, to occur for a brief but finite dura- 
tion in the intermediate time regime, on a CT sinusoidal 
potential, was extended to the case of periodically re- 
versing constant tilts. We have shown the possibility of 
obtaining coherent particle motion interspersed by dis- 
persive motion over many periods of an external ZMSW 
field. The cumulative duration of coherent motion can, 
thus, be extended to a substantial fraction of the total 
journey time, of course, at a cost of making the system 
several times more diffusive. 

MCM acknowledges support from the Abdus Salam 
ICTP, Trieste, where the paper was finalized. 
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